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Abstract 

We study diffraction and interference of indistinguishable particles. 
We consider some examples where the wavefunctions and detection prob- 
abilities can be evaluated in an analytical way. The diffraction pattern of 
a two-particle system shows notorious differences for the cases of bosons, 
fermions and distinguishable particles. In the example of near-field in- 
terferometry, the exchange effects for two-fermion systems lead to the 
existence of planes at which the probability of double detection is null. 
We also discuss the relation with the approach to systems of identical 
particles based on correlation functions. In particular, we shall see that 
these functions reflect, as in noise interferometry, the underlying periodic 
structure of the diffraction grating. 

1 Introduction 

Exchange effects, associated with the peculiar behaviour of indistinguishable 
particles, have been extensively studied in the literature. In recent times these 
studies have followed two principal lines. On the one hand, the experimental 
measurement of correlation functions has corroborated the existence of bunching 
and antibunching effects. On the other hand, some interference phenomena are 
strongly dependent on the distinguishable or undistinguishable nature of the 
particles involved. In this paper we shall be mainly concerned with this second 
line of research. 

In the case of photons, following the seminal work of Hong-Ou-Mandel 
(HOM)[l], the effects of the statistics of the particles have been extensively 
studied for interferences originated by the interaction in a beam splitter. Some 
authors have considered the extension of that approach to the case of massive 
particles. In [2] (and references therein) the behaviour of fermions in Mach- 
Zchndcr interferometers has been analyzed. The proposal presented in [3, 4], 
is also interesting, where an electronic HOM-typc interferometer is used to de- 
tect entanglement. In all the above arrangements the particles can only be 
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detected in the output arms of a beam splitter or interferometer. We analyse 
in this paper whether exchange effects are also present in other interferometric 
arrangements with a continuous spatial distribution. The natural framework to 
discuss this possibility is that of diffraction gratings, which generate continuous 
spatial patterns. 

At this point we must signal the existence of other works by Mandel's group 
where the spatial dependence of two-photon interferences has been analysed 
[5, 6]: two photons generated by down-conversion interfere at a beam splitter 
and detectors are placed in the two output arms. Moving the detectors relative 
to 

the beam splitter, a spatial two-photon interference pattern is observed. At 
variance with our goal, these interferences are generated at a beam splitter and 
not at a diffraction grating. 

As a first approximation to the problem, we consider in this paper some 
simple examples of diffraction and interference by gratings which can be solved 
analytically. This way we can illustrate the main ideas involved in the problem 
without addressing other more technical questions. To be concrete, we shall 
analyse two archetypical systems: 

(i) Diffraction by a single slit of a Gaussian wave packet. This is probably 
the simplest system where the effects can be analysed. We shall evaluate the 
probability detection patterns showing sharp differences between two-particle 
distributions of bosons, fermions and distinguishable particles. 

(ii) Near-field interferometry with periodic gratings. The best known exam- 
ple of this type of arrangement is the Talbot Lau interferometer [7, 8], which 
has been used in a number of interesting studies in matter wave interferometry 
[8, 9]. In this case we shall find a novel characteristic associated with the statis- 
tics of the system. For a pair of indistinguishable fermions, there are planes 
with a null probability of double detection. 

As signaled before, in the other principal line of study of indistinguishable 
systems the focus is on the correlation functions. This line originated in the sem- 
inal work of Hanbury Brown and Twiss (HBT) [10], in which they counted joint 
detections (in two separate detectors) of photons from different chaotic sources. 
The correlations showed that the photons tend to arrive bunched in groups. 
This behavior can be understood within a classical framework if one introduces 
fluctuating phases. In contrast, the opposite trend observed for fermions has 
no classical analogue. Then the only framework able to completely describe 
bunching and antibunching is the quantum one, where these effects can be ex- 
plained in an unified way in terms of constructive and destructive interference. 
In the last years, the correlation functions of free boson [11, 12, 13] and fermions 
[14, 15, 13] have been experimentally obtained, corroborating the existence of 
bunching and antibunching effects. It is natural to ask for the behaviour of 
the correlation functions in interferometry and to see if the presence of inter- 
ference effects can modify the usual picture. A well-suited technique for the 
position-sensitive counting of particles when one is concerned with the spatial 
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dependence of the correlations is noise intcrferometry. This technique converts 
spatial patterns into an interference signal. In particular, in Refs. [12, 15] the 
existence of periodic quantum correlations was shown between density fluctu- 
ations in an expanding atom cloud when released from an optical trap. These 
spatial correlations reflect the underlying ordering in the lattice trap. Noise 
intcrferometry, is useful for identifying quantum phases of ultracold atoms in 
periodic potentials. We shall evaluate the correlation functions in the framework 
(ii) showing that, as in noise intcrferometry, they reflect the periodic structure 
of the underlying arrangement. 

The plan of the paper is as follows. In Sect. 2 we present the basic equations 
and derive some general properties of the interference of two identical particles 
in a diffraction grating. The two analytically solvable models considered in the 
paper are discussed in Sects. 3 (diffraction) and 4 (near-field interferometry) . 
The last one is divided into two subsetions devoted, respectively, to single- and 
multi-mode states. Section 5 deals with the behaviour of correlation functions in 
diffraction gratings. Finally, in Sect. 6, we recapitulate on the principal results 
of the paper and consider the possibility of testing them. 

2 General expressions 

The arrangement we consider here consists in a system of two identical particles 
arriving on a diffraction grating. After the diffraction grating we place detectors 
measuring the interference pattern, which is obtained after many repetitions of 
the experiment. 

We denote by ^>k(x, t) and ip p (y,t) the wavefunctions of the two particles, 
where x and y are the coordinates of the particles and k and p are their wavevec- 
tors or momenta. When the particles are in single-mode states, they refer to 
the wavevectors or momenta of these modes. On the other hand, if we are 
dealing with multi-mode distributions, they represent the mean values of these 
distributions. When the particles arc identical the usual product wavefunction 
tl>] l {'X;t)il) p (y,t), must be replaced by 

-J= (Vk(x, t)i> p (y, t) ± V P (x, *)V>k(y, t)) (1) 

In the double sign expressions the upper one holds for bosons and the lower one 
for fermions. The probabilities associated with this wavefuntion are 

P(x,y,i) = ^ k (x,i)| 2 ^ p (y,t)| 2 + ^ P (x,t)| 2 |^ k (y,t)| 2 

±Re(r p (y, i)V-k(x, i)Vk(y, t)V> P (x, i)) (2) 

We have an interference term that is not present for distinguishable particles, 
according to the standard interpretation of exchange effects as interference ef- 
fects. From now on, in order to avoid any possible confusion by the use of the 
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interference concept in two different ways, we restrict it to the effects induced by 
diffraction gratings. The interference effects associated with the indistinguish- 
able character of the particles will be denoted as exchange effects. 

ft must be remarked that the (anti)symmctrization procedure is only physi- 
cally required when there is a non-negligible overlapping between the two wave- 
functions. If not, the particles must be treated as distinguishable ones. The 
important point is that once the overlapping has taken place, due to the lin- 
earity of the evolution equations, the (anti)symmetric form of the two-particle 
wavefunction persists at subsequent times. 

Several consequences directly emerge from the above expressions: 

(a) For bosons in the same state (p = k), the probability of simultaneous 
detection of two bosons at the same point is P(x, x, t) = 2|?/>p(x, t)\ 4 . The prob- 
ability is proportional to the fourth power of the wavefunction modulus. This 
behaviour was described in models of double ionization by massive particles [16] 
and double absorption of massive particles [17]. At this point a word of cau- 
tion is in order if one tries to experimentally corroborate the above probability 
distribution. The probability in Eq. (2) for different points, x ^ y, can be 
experimentally measured using two detectors placed at x and y. In the case of 
double (boson) detection at the same point x = y, we can only use one detector, 
which must be able to distinguish between single- and double-detection events. 
This type of detector is already available for photons [18, 19], but up to our 
knowledge, not for massive particles. 

(b) For fcrmions the wavefunction of the complete system can have nodal 
points not present for factorizable ones. For example, if the individual wave- 
functions obey (up to a global phase) the relations ^(X, t) = ?/> p (X, t) for a pair 
of points X and Y, these points are nodal ones of the two-particle wavefunction 
(not present in the factorizable case ^(X, t)tp p (Y, t) unless X or Y are them- 
selves nodal points). This property will play a central role in the discussion of 
the detection distributions in near-field arrangements. 

(c) The set of points (X, Y, t) has another interesting property (for bosons 
and fermions). The detection probabilities at these points show a maximal 
deviation with respect to distinguishable ones. A good measure to quantify 
how much the detection rate is enhanced or diminished by the presence of the 
exchange effects is given by the ratio of the detection rates for idistinguishable 
particles and for the same state without exchange effects. It is defined at each 
point and at a given time as 

^(x. y. *) = p 7 7T (3) 

PNE{X,y,t) 

where Pme denotes the probability distribution without interchange effects. In 
the particular case of the set (X, Y, t) we obtain the maximum difference, 

^ (x ' Y ' t} - kMx,t)pkMY,t)P ~ 1±1 (4) 
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which is 2 for bosons and for fcrmions. The rate of simultaneous detections for 
bosons doubles that without exchange effects, whereas that of fermions drops 
to zero. The resemblances with the bunching and antibunching effects are clear, 
although in this case the points X and Y (at a given time t) do not need to 
be very close. The particles were close enough only at a previous time, leading 
to the (anti)symmetrization of the wavefunction, which persists at subsequent 
times. 

In order to go beyond these general properties and to look for other ef- 
fects, we must consider particular arrangements and evaluate the form of the 
wavefunctions in each situation. We do it in the next two sections. 



3 Diffraction 

We consider in this section the diffraction of a pair of indistinguishable particles 
by a slit. The rigorous analysis of the problem leads to the numerical integration 
of Fresnels functions [20, 21]. A simpler approach to the problem is, following 
Feynman [20], to replace the slit by a Gaussian slit. With this approximation 
we have an analytically solvable problem. The wave packet generated by this 
slit with soft edges has also a Gaussian profile in the direction parallel to the slit. 
The movement in the perpendicular direction is assumed to be unaffected. As 
usual in the treatment of the problem, we consider a two-dimensional problem, 
neglecting the vertical axis. Moreover, as the movement in the perpendicular 
axis is almost unaffected, the relevant physics is contained in the parallel axis 
and the problem reduces to an one-dimensional one. 

The simplest way to study multimode states is to assume that the wavevec- 
tor distribution of each particle is a Gaussian one [21]. Using the notation of Ref. 
[21] the one-dimensional Gaussian distribution is f(k) = (Air) 1 / 4 a -1 / 2 cxp(— (fc— 
fco) 2 /2cr 2 ) with a the width of the distribution and ko its central value. The 
wavefunction is obtained by superposing a set of planes waves with that distri- 
bution ((271-)- 1 / dkf(k)exp(i(kx - E k t/h)) with E k = h 2 k 2 /2m) [21]: 

,, ^ „,r. ( -o 2 {x - vt) 2 + ik Q (2x - vt) + iha A x 2 t/m\ ... 
ip(x, t) = C(t) exp i '- K — '- '- — (5) 



with v = hko/m, 



and 



/*(*) 



C(i)= 7 r- 1 / 4 ( 1 + — (6) 
(j ml 



1 ' ^ ) 



Note the absence with respect to the formula in [21] of a multiplicative coefficient 
2 in the last term of the numerator in the exponential. 
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As required by our assumption of a Gaussian slit, the wavefunction has the 
form (up to some phase terms) of a spatial Gaussian packet with time-dependent 
width a' 2 + (h 2 t 2 a 2 /m 2 ). 

Using the expression of the wavefunction, it is immediate to evaluate the 
probability distribution for two indistinguishable particles. We assume both par- 
ticles to be in states of the type (5), with the same width and central wavevectors 
fco and pq. The final result is 

P(x, y, t)/\C(t)\ 4 = l -[e A + e B ± 2e^ B )l 2 cos (J^{x - y)(k - p )) ) (8) 

with 

2a 2 

A = -^{{x-v k tf + (y~ V) 2 ) (9) 

and B given by a similar expression with the interchange of v k and v p . 

From the above expressions it is clear that the detection probabilities for 
a pair of indistinguishable particles show notorious differences with respect to 
those of distinguishable ones. We illustrate this point numerically in Fig. 1, 
where we compare the probability of simultaneous detection for bosons, fermions 
and distinguishable particles. In order to simplify the presentation, we fix one 
of the points, x — x , and evaluate the detection probability for arbitrary y. 
We choose for xo, at a fixed time t, the value xo — v k t which maximizes the 
detection probability in the first exponential. The probability in the case of 
distinguishable particles is Pdis{x, y, i)/|C(i)| 4 = e A (x = x ) = cxp(— 2<r 2 (y — 
v p t) 2 / n{t)). Taking M/m = 0.1, k = 1 = —po and a width of the distribution 
a 2 = 0.125 we have the following distribution: 
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Figure 1: In the vertical axis we represent, in arbitrary units, p — 
P(xo,y,t)/\C(t)\ 4 . In the horizontal axis we have the coordinate y, also in 
arbitrary units, and the point corresponding to y = v p t. The red, black and 
blue curves correspond, respectively, to bosons, distinguishable particles and 
fermions. 

Several consequences emerge directly from the figure. As expected, the dis- 
tribution in the distinguishable case retains a Gaussian form centred around 
the point y — v p t. In contrast, for both fermions and bosons, the mathematical 
form of the curves departs from the Gaussian profile showing the typical form of 
an interference pattern. The maximum probability for bosons is approximately 
at the same position of the peak of the Gaussian for distinguishable particles 
and doubles its value. At the same position we have a null-probability point for 
fermions. 

All the above development has been done in terms of the diffraction by a 
slit. However, these expressions also describe two particles that interact (there 
must be a strong overlapping between them) at a given moment and then con- 
tinue a free evolution. We can measure the positions of the particles after the 
interaction, which are given by Eq. (8). Note that in this case that expression 
is exact, whereas for the description of the diffraction by the slit, it is only an 
approximation. 

4 Near-field interferometry 

After considering an example of diffraction, we now move to a purely inter- 
ferometric arrangement, which illustrates other types of effects associated with 
the interchange terms. In particular, we focus on near-field interference in the 
eikonal approximation [8] , which can be tackled analytically (for the limitations 
of the eikonal approximation see [22]). One of the best known examples of 
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near-field arrangement is the Talbot Lau one [7], which has important applica- 
tions in matter wave interferometry [8, 9]. In order to emphasize on the main 
physical ideas of the problem without excessive mathematical technicalities, in 
a first step we shall only consider particles in single-mode states, showing later 
that the main result persists (under adequate restrictions) in the more realistic 
framework of multi-mode states. 

4.1 Single- mode states 

We consider a plane wave ip a — exp(ifcz) which incides on a grating placed in 
the plane perpendicular to the z-axis. If the grating function is a periodic one, 
the transmission function is given by 

T(x) = ^2 A n exp(2mnx/d) (10) 

n 

with n an integer, d the grating period and A n = \A n \ cxp(i£„) the coefficient 
of the term n. The wavefunction at a distance L behind the grating is [8]: 




with Afc the wavelength of the particle and x the coordinate of the axis in which 
we measure the interferences. 

We consider now the normalization of this wavefunction. From the expres- 
sion 

\^{X)\ 2 = E A n A m ™P [ ~ ^ ~ j ™P ( — ~ J (12) 

we see that \ip k (x)\ 2 dx is not bound. In effect, the integration of the 
constant terms AqAq, A\A\, ■ ■ ■ docs not give a finite value. One possibil- 
ity to overcome this difficulty is to do the normalization in a relative sense, 
Af- 1 f™ x \^(x)\ 2 dx = 1, with TV = \A n \ 2 ) dx. An alternative form of 

normalization is to consider the integration over the period d, 4 Jq dx\ip k (x)\ 2 = 
1. Taking into account that dx ex.p(2irinx / d) — d5 n o, we have ^ dx\ip%(x)\ 2 — 
J2 n \A n \ 2 = A 2 , and the normalization condition is obtained dividing by A. It 
is simple to see that the same condition is reached using the relative normal- 
ization condition. The normalized wavefunction after the diffraction arrange- 
ment is given by ip%(x)/A, or equivalently by Eq. (11) with the replacement 
A t -> a, = Aij A (a t = \A t /A\ex P (i£ n )). 

As anticipated in Sect. 2 let us consider the points X where the condition 
ip%(X) = ipp(X) holds (up to a global phase). In this example, the problem 
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is stationary and we do not need to include explicitly the time variable. We 
introduce the new parameter 

ttL 

<t>k P = "^"(Afe - A p ) (13) 

For all the initial wavevectors fulfilling the condition cpkp — 27m*, with n* an 
integer, we have that at all the points X contained in the plane defined by 
z = L, the wavefunctions obey the relation 

^(X) = e^ L ^(X) (14) 

From now on, we restrict our considerations to fermions. The antisymmetrized 
wavefunction at that plane becomes null: 

-L(^(X)^(Y) ^(X)^(Y)) -+ 

-Le^ L (^(X)^(Y) - tf{X)tf(Y)) = (15) 

For all the initial wavevectors for which the relation 4>k p = 2im* holds, we have 
that in the plane z = L, there cannot be two-fcrmion detections. Conversely, 
once the values of Afc and X p are fixed, there are always planes in which the 
phenomenon can be observed. These planes are given by z — L n , with L n = 
2nd 2 / (Afc — A p ) and n any integer. All the points in these planes are nodal points 
of the two-fermion wavefunction. 

Note that this result does not forbid one-fermion detections at these planes. 
In effect, the one-fermion detection probabilities are given by the reduced de- 
tection probabilities, V(X) — J P(X,y)dy. The integration on y contains any 
y placed on any plane, not only the Y's contained in the nodal plains. The 
contributions of the points in the nodal plains are null, P(X, Y) = 0, but the 
contributions P(X, y) with y outside these special planes can be different from 
zero. The sum of all these contributions can be, in principle, different from zero 
leading to a non-null one-fermion distribution at the nodal planes. 



4.2 Multi-mode states 

We move now to the more realistic case of non-monochromatic states. We show 
that the result of the previous subsection is still valid for narrow enough initial 
wavepackets. 

If the momentum distribution of the wavepacket is f(k), the wavefunction 
at a distance L behind the grating becomes 

^ {x ) = ^A n F n exp(^) (16) 
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with 



F n = J dkf(k)e lkL e { - l7Tn2LXk/d2) ~ 

e ik L e (~i7rn 2 L\ ka /d 2 ) J ^g-( fe - fe o) 2 /2cr 2 e %{k-k a )L e (-iTtn 2 L(k -k) /d 2 kk ) ^) 

where we have adopted for the mode distribution a Gaussian curve with the 
central wavevector fco- In order to simplify the notation, we have dropped the 
coefficients associated with its normalization. 

The two imaginary exponentials inside the integral have a joint argument 
L(k — fc )(l + (irn 2 /d 2 kk )). When the condition d 2 kk a » im 2 holds, the 
integral simplifies to I(k ) = J dke~( k ~ k °^ /2<T e *( fe - fe o)£_ L e t us consider the 
circumstances under which this approximation can be justified. Three condi- 
tions must be fulfilled. (1) n must be small. For large values of n the right-hand 
side of the inequality will be very large and the term containing it cannot be 
neglected. Thus, the summation in (16) must be truncated at some value no. 
This truncation is plainly justified because most of the behaviour of the system 
can be described taking into account only the lower terms in the summation. As 
a matter of fact, in the representative experiment in [8] the interference signal 
is essentially determined by the and ±1 components only. (2) For values of k 
close to the central one, ho, we must have d 2 kko ~ d 2 k^ » Trn 2 ,. This condition 
can be fulfilled by restricting our considerations to this range of central wavevec- 
tors. (3) As a consequence of (2), the distribution must be sharp enough around 
fco in order the contribution of the elements with \k\ -C |fc | to be negligible. 

With this approximation Eq. (16) becomes: 

^k ( x )= J2 A n I{k y koL e { - m2LXk « /d2) e {2mx/d) (18) 
M<l™ol 

Let us consider now another distribution with the same Gaussian form (the same 
width), but centered around another value po- Then /(fco) = I(po), and we are 
at the same position of the previous subsection: when the relation <j)k oPo = 
27rn* holds we have that in the plane z = L the wavefunctions obey Vfe (A) = 
e i(fc -Po)i^£ (x). 

Within the range of validity of this approximation, the result obtained in 
the previous subsection for a single-mode initial wavefunction can be extended 
to its multi-mode counterpart. 

5 Correlation functions 

As signaled in the introduction, another fundamental line of research in identical 
particles is based on correlation functions. We analyse in this section whether 
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the presence of the interference device modifies the behaviour of the functions 
found in free systems or optical lattices. 
The correlation function is defined as 



1 f d 

C(rj) = - / P(x, x + r))dx 
" Jo 



(19) 



To calculate this expression we must evaluate the two-particle probability dis- 
tribution. In order to simplify the calculations, we restrict our considerations to 
the case of particles described by single-mode plane waves. As we have seen in 
the previous section, the result obtained for a mono-chromatic particle can be 
extended to sharp enough wave packets. We analyse the correlations in a fixed 
plane z = L. The two detectors measuring the correlations must be placed in 
that plane. Then it is not necessary to include the L parameter explicitly as a 
superscript of the wavefunction. 

The modulus of the initial state is: 

„ /27t(to — n)x TrL(m 2 — n 2 )\ k 

\Mx)? = 1 + 2 22 Ka m \ COS ( — ^ " ^2 — +Zm-Zn 

n,m ^ 

(20) 

In order to complete the expression for the probability, we must evaluate the 
crossed or interference term 

a*n a *m a ras X 

2m[{s - n)x + (r - m)y] \ f iwL[X p (s 2 - to 2 ) + A fe (r 2 - n 2 ) 



exp exp 



d 1 <\ d 2 

E \a n a m a r a s \ x (21) 



( 2tt[{s - n)x + (r - m)y] _ nL[X p (s 2 - to 2 ) + A fc (r 2 - n 2 )} | 

C0S \ d d 2 +«r,w 

where £, r ,s,n,m = £r + £s — £n — £m ■ Adding the direct and crossed terms, we 
obtain for the probability 

P(x,y) = l + 

m>n 

E K a ™K C0S< £m,n( x ) +cos^„(x) +cos^„(y) +cos^„„(y)) + 

m,n 
7ii>n,r>s 

2 E \a n a m a r a s \(cos ip^Jx.) cos ipP s (y) + cos^ „(x) cos^ s (y)) ± 

E |a n a m a r a s | x (22) 
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/ 2tt[(s - n)x + (r - m)y] nL[\ p (s 2 - m 2 ) + \ k (r 2 - n 2 )] \ 

° 0S { d d 2 +Sr,s,n,mJ 

where <p k m Jx) = 2 -I^ph _ *Hm>- n >)x k + ^ _ ^ Now ^ the correlation 
function can be obtained from Eq. (22). The calculation is simple but lengthy. 
There are four types of contributions: (i) constant terms, (ii) terms containing 
cos¥>£j n (x), (hi) terms with cos ^J, (a;) cos <p% r (y) and (iv) those contained in 
the crossed term. Let us evaluate them separately. 

(i) For constant terms we have, g J ctedx — cte. 

(ii) By direct calculation we have, g cos ip k nn (x)dx — = ^ J Q d cos ip^„(i+ 
rj)dx. The terms with this form do not contribute to the correlation function. 

(iii) Using standard trigonometric relations and integrals we have that the 
integral ^ cos iPm,n( x ) cos Vr.s^ + v)dx equals when m-n^i(r-s) and 



2 COS [ " d ' + &sr) + 2 C0S { d + (23) 

otherwise, with (f>% nsr = -$[(m 2 - n 2 )\ k + (r 2 - s 2 )\ p ] + £ m - - £ s + 

(iv) For the crossed term, we must distinguish when s — n + r — m is zero or 
different from zero. When it is not zero, we have again the situation (ii) and its 
contribution is null. On the other hand, when s — n + r — m = 0, the integral 
of the cosine gives 

(^^ + e*„) (24) 

with e k n P mrs = -^[(r 2 - n 2 )\ k + (s 2 - m 2 )\ p ] + &, s , n , m . 

Combining all these expressions, we obtain the correlation function: 

m>n,s>r 

C(rj) = 1 + 2 ^ \a m a n a s ar\ x 

n,m,r,s 



COS 



V d ' T mnsr I 2 V d ' > nnsr 

s — n+r— m— 



j>2 \a m a n a s a r \ cos f 2 ^ - + 0^„^ (25) 



n.m,r.s 



The correlation function depends on the parameters of the grating, dj and (i 
and on the wavevectors of the incident particles (through <fim nsr and 0„^ rs 
both functions of Afe and X p ). This behaviour is to be compared with that of 
free particles, for which the correlation function depends on the wavevectors in 
the form Cf ree (rf) — 1 ± cos((p — k)rj) (for plane waves exp(ikx) and exp(ipy)). 
The spatial periodicity of C/ ree is determined by p — k. For interfering particles, 
however, the periodicity is only ruled by d/n (with n any integer), a function 
of the parameter of the grating d. In the correlation function we have all the 
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periods (d/n) contained in the diffraction grating. The wavevectors only deter- 
mine (in conjunction with d and the coefficients £i) a phase for each term. We 
conclude that, as in noise interferometry, the correlation functions reflect the 
underlying structure of the periodic grating. 

A particularly clear illustration of the behaviour of the correlation functions 
is obtained in the particular case that only oti and a are different from zero. 
Moreover, for the sake of simplicity, we assume the two coefficients to be real. 
The correlation function in this case is (using a\ + a 2 , = 1) 

C(r]) = (1 ± 1) ± 2a 2 a 2 (-l + cos^ fep ) + 2a 2 a 2 (±l + cos</> fcp ) cos {^j^ (26) 

Now there is only one periodicity, d. The wavevectors only determine (in con- 
junction with a and a\) the amplitude of the oscillations. Moreover, the de- 
pendence is in the form \ — |, instead of p — k. 

For short distances, r\ -> 0, we have C F {r\) — %-(^ L ) 2 ?7 2 (Co = ^a\a\{\ - 
cos^fe p ))and Cb(v) — Cb(0) <~ — aQa 2 (^ L ) 2 (l + cos (/^p)^ 2 that is the usual 
dependence on r? 2 , but with different coefficients. 



6 Discussion 

In this paper we have analysed the behaviour of exchange effects in diffraction 
and interference. We have found the existence of some distinctive characteristics 
of this type of interferometry: modifications of the probability distributions 
with respect to those of distinguishable particles, strong increase or decrease 
of the double detection rates for some pairs of points (similar to bunching or 
antibunching, but without need of closeness between the particles) and existence 
of planes where the double detection of fcrmions is forbidden. On the other 
hand, the correlation functions show a behaviour similar to that found in noise 
interferometry, reflecting the periodic structure of the diffraction grating. 

Now briefly address the question of how the above distinctive characteristics 
of continuous interferometry of identical particles can be tested experimentally 
As signaled before, the interchange effects are present (and persist in the sub- 
sequent evolution of the system) when there is a non-negligible overlapping 
between the wavefunctions of the two particles at a given time. Then the key 
step is to obtain a non-negligible overlapping of the two particles at the time 
they reach the diffraction grating. This is equivalent to a careful preparation 
of the times of arrival of the two particles. In other words, the peaks of the 
two distributions must reach the grating at the same time. Because of recent 
progress in matter wave interferometry, interferometry of heavy molecules and 
the existence of massive systems with high degrees of coherence, these effects 
seem to be accessible to experimental scrutiny. Conversely, the presence of the 
effects described in this paper could be used as a test of the closeness of two 
particles under some specific preparation. 
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It must also be noted that because of the interaction of the beam with the 
grating, there can be absorption and backscatter processes which result in cases 
with zero or one particles arriving at the detectors. We can overcome this 
difficulty by considering a postselection process in which only events with two 
detections are taken into account. 

In this paper we have restricted our considerations to two simple examples 
that can be treated analytically. Other types of diffraction gratings must be 
studied in order to see if (the same or different) exchange effects are also present. 
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